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ABSTRACT 

We derive an analytic expression for the mean magnification due to strong 
gravitational lensing, using a simple lens model, a singular isothermal sphere 
embedded in an external shear field. We compute separate expressions for 2- 
image and 4-image lensing. For 4-image lensing, the mean magnification takes 
a particularly simple form: (/i 4 ) = ^tj~2) , where 7 is the external shear. We 
compare our analytic results to a numerical evaluation of the full magnification 
distribution. The results can be used to understand the magnification bias that 
favors the discovery of four-image systems over two-image systems in flux-limited 
lens surveys. 



1. Introduction 

There have been several systematic surveys for examples of multiple-image gravitational 
lensing of extragalactic sources, both at optical wavelengths (e.g., Maoz et al. 1993; Surdej et 
al. 1993; Gregg et al. 2001) and radio wavelengths (Burke et al. 1993; King et al. 1999; Winn 
et al. 2000; Browne et al. 2002). The interpretation of the number of lensed sources that 
are detected in a flux-limited survey depends critically on the factors by which the sources 
are magnified. Lensed objects are drawn from a parent population with smaller fluxes than 
their unlensed counterparts. One therefore samples a fainter portion of the source luminosity 
function and a more distant range of source redshifts. This well-known "magnification bias" 
also affects the relative frequencies of different lens morphologies. Four-image lenses are 
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more highly magnified than two-image lenses, and are subject to a stronger magnification 
bias. 

Early studies of lens statistics used lens models with circular symmetry, such as singular 
isothermal spheres, for which the magnification bias is easy to calculate (e.g., Turner, Os- 
triker, & Gott 1984). However, symmetric models only produce 2-image lenses and cannot 
be used to interpret the statistics of each morphology separately. More recently, numerical 
studies have been performed with more realistic elliptical lens models (e.g., King & Browne 
1996; Wallington & Narayan 1993; Keeton, Kochanek, & Seljak 1997; Rusin & Tegmark 
2001). While numerical calculations are certainly valuable, they tend to obscure the roles 
of competing contributions to what is being calculated. Analytic forms for such things as 
luminosity functions and magnifications may be no more accurate, but they produce results 
that are more readily understood. 

With this motivation we examined a number of simple elliptical lens models and found 
one for which the mean magnification can be computed analytically, for both 2-image and 4- 
image lensing. The model is a singular isothermal sphere in an external shear field (SIS+XS), 
which is physically motivated and widely used. 

In the next section we briefly review the mathematics and the qualitative features of 
lensing by non-circular potentials. In § 3, we specialize to the SIS+XS and derive the analytic 
expressions for all of the relevant quantities, including the mean magnification. Finally, in 
§ 4 we compare our result to the full numerical evaluation of the magnification distribution, 
and remark on applications to the interpretation of lensing statistics. 



Lensing theory is well developed and has been expounded in detail by Schneider, Ehlers, 
& Falco (1992), among others. The particular problem of lensing by elongated potentials has 
been lucidly described by Blandford & Narayan (1986). Here we briefly review only those 
concepts that are required to follow the calculations presented in § 3. 

We use the conventional definition of the lens potential, 



in which $ is the Newtonian gravitational potential of the galaxy; r is the sky-plane angular 
displacement; z is the line-of-sight coordinate; and Di, D s , and Di s are the angular-diameter 
distances to the lens, to the source, and between lens and source, respectively. With this 



2. Review of lensing by elliptical potentials 
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definition, the correspondence between source position s and image position r is 



s = r — Vr^(r). 



(2) 



The inverse magnification of a particular image % is given by 
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which is negative when the image and source have opposite parity. For a survey in which 
sources are selected by total flux, the relevant quantity is the total magnification /^(s) = 
^2 (s) | , the sum of the unsigned magnifications of all images of a source located at s. 

Most of the observed galaxy lenses can be described (at least qualitatively) using cen- 
trally concentrated, singular potentials. If the isopotential contours are only moderately 
elliptical, such potentials produce either 1, 2, or 4 images of a background source, in loca- 
tions that can be understood with reference to the curves plotted in the top two panels of 
Figure 1. The curves are plotted for the SIS+XS, but the following discussion applies to 
many singular elliptical potentials. 5 

For a source outside the radial caustic 6 , one image is produced outside the 1-2 transition 
locus. A source on the radial caustic corresponds to one image on the 1-2 transition locus 
and another image at the origin, with zero flux. If the source moves from the radial caustic 
to the astroid caustic, the outer image moves to the outer 2-4 transition locus, while the 
inner image acquires a nonzero flux and moves to the inner 2-4 transition locus. On the 
astroid caustic, three images are produced: one on each of the 2-4 transition loci; and one 
on the critical curve, with infinite flux. Inside the astroid caustic, the latter image splits into 
two images with finite flux, and the result is one image in each of the regions marked "4" in 
Figure 1. The caustics and critical curve are well known to those who study gravitational 
lensing, but the transition loci receive less attention than they deserve — in particular, they 
are crucial for our calculation of mean magnification. 

For more elongated potentials, two cusps of the astroid caustic pierce the radial caustic 
and are often referred to as "naked cusps." The relevant curves in the source plane and 
the image plane are shown in the bottom panels of Figure 1. A source inside a naked cusp 



5 Non-singular lens potentials produce an odd number of images. In such cases, the image plane has an 
additional critical curve and transition locus surrounding the origin. 

6 Strictly speaking, this curve is not a caustic, because the potential is singular. For this reason, the curve 
has also been designated the radial pseudo-caustic (Evans & Wilkinson 1998) or cut (Kormann, Schneider, 
& Bartelmann 1994). 
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Fig. 1. — Top. Source plane (left) and image plane (right) for the SIS+XS with b = 1 
and 7 = 0.2. The curves discussed in § 2 are labeled. The regions are labeled by image 
multiplicity, e.g., a point source in a region of the source plane marked "4" corresponds to 
four points in each of the four regions of the image plane marked "4." Bottom. Same, but 
for 7 = 0.45. 
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corresponds to three images located on the same side of the origin as the source — one in 
each of the three regions marked "3" in Figure 1. 

If the source position is random, the cross section for producing 4 images, <t 4 , is the total 
area within the 4-image region of the source plane (marked "4" in Figure 1). The 2-image 
cross-section is the area of the 2-image region, and the 3-image cross-section is the area of 
the 3-image region. The cross-section for multiple imaging is a m = a 2 + <r 3 + <t 4 . 

The mean magnification for multiple-image lensing is the average value of the total 
magnification for all source positions that result in more than one image: 



where the integral of dr extends over the region in the image plane in which multiple images 
are found. Eq. 4 holds because the Jacobian for the s — > r transformation is equal to the 
magnification. Alternatively, one can visualize the calculation as follows. Because lensing 
conserves surface brightness, a source consisting of uniform brightness in the region formed 
by the union of the caustics will result in an image of uniform brightness within the union 
of the transition loci. The mean magnification for multiple-image lensing, which we denote 
(fj, m ), is the ratio of the image area to the source area. The mean magnification of 4-image 
lenses, denoted (/X4), is the ratio of 4-image region of image space and the 4-image region of 
source space. The mean magnification of 2-image and 3-image lenses can be visualized in a 
similar manner. 

Of the approximately 100 observed cases of multiple-image gravitational lensing by 
galaxies, almost none require naked cusps for explanation. The one possible exception is 
APM 08279+5255 (Lewis et al. 2002). Because naked-cusp systems are apparently rare in 
nature, and because the calculations described below cannot be done analytically for naked- 
cusp systems, in this paper most of our attention will be on systems in which the astroid 
caustic is interior to the radial caustic. 

In that case, to find the cross-sections and mean magnifications for 2-image and 4-image 
lensing, it is sufficient to calculate the areas enclosed by the radial caustic (A r ), astroid caustic 
(A a ), the 1-2 transition locus (At), and the inner and outer 2-4 transition loci (Ai and A Q ). 
As long as the astroid caustic is contained entirely within the radial caustic, we have: 
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3. Analytic expressions for the SIS+XS 

We concentrate on the case of a singular isothermal sphere in an external shear field 
(SIS+XS), for which the potential is 

2 

ip(r,6) = br + - 1 — cos 20. (6) 

Here (r, 9) are polar coordinates in the image plane, b is the Einstein ring radius (which 
sets the mass scale of the lens), and 7 is the magnitude of the external shear. The SIS is an 
idealization of a spherical dark-matter halo with a flat rotation curve, and the external shear 
is a non-circular perturbation. This potential can be thought of as a truncated multipole 
expansion of the lens potential (see, e.g., Kovner 1987; Schneider & Weiss 1991; Trotter, 
Winn, & Hewitt 2000). The SIS+XS has been used extensively in the literature to model 
observed lenses, often as a starting point for more complicated models (see, e.g., Kochanek 
1991; Keeton, Kochanek, & Seljak 1997; Lehar 1997; Schechter et al. 1997). 

By computing the second derivatives of this potential and using Eq. 3, we obtain the 
magnification 

/i" 1 = 1 - 7 2 - h - (1 - 7 cos2#) . (7) 

We now derive expressions for all the relevant curves and the areas enclosed by them. 
We identify the radial caustic (x r ,y r ) by requiring that it map to the origin under Eq. 2, 
obtaining 

x r = —b cos 9 

y r = — &sin#, (8) 

which is a circle of radius b, enclosing an area A r = nb 2 . 

Next, we identify the 1-2 transition locus by reverse- mapping the radial caustic back 
into the image plane, and finding the non-zero solution r t (9 t ): 

x r = —b cos 9 — (1 — 7) r t cos 9 t — b cos 9 t 

y r = — b sin 9 = (1 + 7) r t sin 9 t — b sin 9 t . (9) 
We eliminate the parameter 9 by squaring and adding the equations, leading to the solution 

«*>-»( iXr£i> )- <io) 
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The area enclosed by the 1-2 transition locus can be expressed in terms of definite integrals 
that are analytically tractable: 
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To identify the critical curve r c (6 l c ), we require the inverse magnification ji 1 given by 
Eq. 7 to be zero, obtaining 

&(1 - 7 cos2fl c ) 

= • ( 12 ) 

We obtain parametric equations for the astroid caustic (x a ,y a ) by requiring the astroid 
caustic to map to the critical curve under the lens mapping, Eq. 2: 
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The maximum value of x a is x* = and the area within the astroid caustic is 



y a dx a 
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We note that the astroid caustic is completely interior to the radial caustic for 7 < 



Under the lens mapping, the astroid caustic has 2 images apart from the critical curve: 
the inner 2-4 transition locus, r^Oi), and the outer 2-4 transition locus, r o (0 o ). In order to 
solve for these transition loci, which we refer to generically as r 4 (^ 4 ), we insert the parametric 
equations 13 into Eq. 2: 



COS Or 

1 + 7 
267 



sin 3 Or 



7 



(1 — 7) r 4 cos 4 — b cos #4 
(1 + 7) r 4 sin O4 — b sin O4. 



(15) 



After eliminating r 4 from this pair of equations, and defining t = cos C , we obtain 

(cos# 4 -t) 2 [cos 2 # 4 + 2t(l-t 2 )cos# 4 -t 4 ] =0, (16) 



- 8- 



which has three solutions: the double root cos 9^ = t corresponding to the critical curve, and 
the 2-4 transition loci 

cos ijO = t (t 2 - 1 ± Vt 4 - t 2 + l) . (17) 

We insert these solutions for cos #4 into the first expression of Eqs. 15 to obtain parametric 
equations for the 2-4 transition loci: 

2 7 t 3 



r 4 (t) 



1-7 



1 - 



(18) 



(1 + 7) cos e A (t)_ 

It is possible to find the area within each locus separately — but for the calculation of mean 
magnification, the key quantity is the area between the 2-4 transition loci, which is (after 
lengthy algebra) 

A o -A i = 4x-J (rl - r\)d9 = ^(Jj + J 27 ), (19) 
where I\ and 1% are defined as 



7T / n COS6 1 



/, = i 

7T 



cos COS 2 



cZ0, (20) 



and ti t0 (6) are the inverse relations of Eq. 17. It is straightforward to prove 7 2 = by 
demonstrating that the integrand is antisymmetric about 9 = j (see Finch 2001), and I\ is 
a dimensionless number independent of b and 7. Numerically we find I\ = 1.35111, leading 
to the result 

47r6 2 J l7 176 2 7 



- A - 7; -57J ~ t; -0X2- ( 21 ) 



Now that all of the relevant curves and areas have been calculated, we compute cross 
sections and mean magnifications using Eqs. 5. We emphasize that these results are only 
valid for 7 < |, when the astroid caustic is completely interior to the radial caustic. To find 
mean magnifications for 7 > |, one must determine the intersection points of the various 
curves, for which we do not have analytic solutions. In § 4 we present numerical results for 
this case. 

The SIS+XS is the only example we have found that has all the qualitative features 
desired of a model for an elliptical potential and for which the areas bounded by the transition 
loci are all analytic. An example of another commonly-used potential for which we do not 
have analytic formulas for all the areas is 

i/j(r,9) = br + 7&rcos2#, (22) 
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which can be thought of as a lowest-order multipole expansion of the singular isothermal 
ellipsoid (SIE) density distribution (see, e.g., Kormann, Schneider, & Bartelmann 1994). 
The astroid caustic is contained completely within the radial caustic for 7 < |. For this 
potential we found analytic expressions for all the curves except the 1-2 transition locus, 
using a procedure similar to the one described above. We were therefore unable to compute 
(/i 2 ), but all the other quantities are listed in Table 1. 



4. Discussion 

To put the analytic results in context, we numerically computed the full magnification 
distribution for multiple imaging by the SIS+XS, denoted P m (/j,). We parameterized the 
points within the radial caustic and computed the total magnification of each point, by 
inverting Eq. 2 to find the images and Eq. 7 to derive the magnifications. We performed this 
procedure separately for the 2-image and 4-image regions of the source plane, so that the 
magnification distribution for each morphology (P 2 and P 4 ) could be computed separately. 
The results are plotted in Fig. 2 for the case 7 = 0.2. 

The general features of the distribution can be readily understood. The minimum 
magnification occurs when the source is on the x-axis of the radial caustic, giving 

2 

^,min = (1 + 37)(1 _ 7) - (23) 

There is a discontinuity at the minimum magnification for 4-image systems, which occurs 
for a source at the origin: 

2 

A*4,min = 7\ oT- (~4) 
7 (l-7 2 ) 



Table 1. Cross sections and mean magnifications 



i>(r,9) 7max cr 2 <t 4 (^2) <A*4> 



i -^j(i-v)-i 17 8Ji 

1_5£) ( 1_ 7 2) 3 7 (l-7 2 ) 

br + jbr cos 29 \ tt6 2 (l - f 7 2 ) 6vr6 2 7 2 ••• ^ 



*The number Ji is defined by Eq. 20, and has the approximate value 1.35111. 
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Fig. 2. — Probability distribution of magnification by the SIS+XS model, shown here for the 
case 7 = 0.2. The 2-image distribution (P 2 ) and the 4-image distribution (P 4 ) are plotted 
separately. The minimum magnifications, given by Eqs. 23 and 24, are marked. 



The high magnification tail is produced by sources in the vicinity of the astroid caustic. It 
can be shown generically that the asymptotic probability for caustic-dominated events varies 
as yU~ 3 (see, e.g., Blandford & Narayan 1987; or Schneider, Ehlers, & Falco 1992, p. 319), 
and indeed the logarithmic slope of P4 as computed numerically is —3.1. The downward 
kink between /i = // 2 ,min an d A^min occurs at the minimum magnification for a source just 
outside the astroid caustic. 

The number of multiple-image lenses in a flux-limited sample can be determined from 
the magnification distribution. As described in detail by Turner, Ostriker, & Gott (1984), in 
a survey of objects with intrinsic luminosity function <f>(L, z) and a flux limit corresponding 
to an intrinsic luminosity L\i m , the bias factor B by which multiple-image lenses at a given 
redshift z are over-represented is 

J Lmin 0(L, z)dL 

The bias factors for 2-image and 4-image lensing can also be computed separately, by replac- 
ing P rn with P 2 or P4. For the common case that the luminosity function is approximated 
by a power law, 

^=^Am)"' (26) 

the bias factor reduces to B — ^ / u /3_1 ), which is simply the mean magnification for (3 = 2. As 
pointed out by Rusin & Tegmark (2001), the value (3 = 2 happens to be a good approximation 
for the radio source populations that have been searched in recent lens surveys (JVAS, King 
et al. 1999; CLASS, Browne et al. 2002; Winn et al. 2000). Our analytic results are therefore 
directly relevant to this case, in which previous calculations were performed by numerical 
integration and Monte Carlo simulation. 

For example, the JVAS and CLASS surveys have found roughly equal numbers of 2- 
image and 4-image lenses, a result that various authors have struggled to interpret because 
it seems to require a magnification bias for 4-image lensing that is larger than expected 
(King & Browne 1996; Keeton, Kochanek, & Seljak 1997; Rusin & Tegmark 2001). A 
resolution considered by these authors is that the matter distribution in the lens galaxies is 
systematically more elongated than the observed light distribution of elliptical galaxies. 



We can estimate of the required value of 7 from our calculations of mean magnification. 
The ratio of 4-image lenses to 2-image lenses, including the effect of magnification bias, is 
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For the SIS+XS model and a (3 = 2 luminosity function, the result is 

iV 4 ct 4 (/i 4 ) ii7 



iV 2 a 2 (/i 2 ) (i_ 7 2)_/ l7 ' 



(28) 



This result is rather specialized, depending as it does on the assumption of a particular lens 
model and source luminosity function, but it has the virtue of being analytic. To obtain 
equal numbers of 4-image and 2-image systems we require 7 = 0.32, for which 02/04 and 
B4/B2 are both approximately 5. 

For 7 > I one must also consider the possibility of 3-image systems produced by naked 
cusps. In this case the computation of areas in the image plane must be done numerically. 
Figure 3 shows the results of a numerical calculation of the relative areas of the 2-image, 
3-image, and 4-image regions of the image plane, for < 7 < 1. As before, if f3 — 2 then 
this plot can also be interpreted as the expected fractions of 2-image, 3-image, and 4-image 
lenses discovered in a flux-limited survey. For 7 < | the results are identical to our analytic 
formula. For 7 = |, naked-cusp images are possible but occur with vanishing probability; 
they constitute 10% of the lenses at 7 0.45 and they dominate the sample for 7 > 0.61. 

In more generic settings, the full magnification distribution is required. The 4-image 
magnification distribution, at least, can be approximated by -P 4 (/x) oc /i~ 3 for /i > /i 4m i n . 
The mean of this approximate distribution is 

J . ft 3 d/i 1(1 -T) 

p 4,mm 

which has the same dependence on 7 as the exact expression (see Table 1) and differs only 
by a constant fraction of 4/3.6 or 11%. 



5. Conclusion 

We have described a geometrical method for computing the mean magnification for 
multiple imaging by non-circular lens potentials. For the special case of a singular isothermal 
sphere embedded in an external shear field, we have derived analytic expressions for the mean 
magnification for 2-image lensing and 4-image lensing. The results can be used, in certain 
circumstances, to compute the magnification bias for 2-image and 4-image lensing. 



We thank Nick Morgan, Jackie Hewitt, and David Rusin for helpful discussions. We 
gratefully acknowledge the support of the Graduate Student Office at MIT (T.K.F.), the 
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Fig. 3. — Numerical evaluation of the relative areas of the 2-image, 3-image, and 4-image 
regions of the image plane, for the SIS+XS model. Under the assumptions described in § 4 
these curves also give the ratios of 2-image, 3-image, and 4-image systems expected in a 
flux-limited survey. 
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